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INTRODUCTION 


This  annual  summary  report  describes  research  performed  from  1  August  1991 
through  14  August  1992  with  support  from  ONR  under  grant  N00014-89-J-1003. 
The  following  projects  are  discussed  in  this  report: 

I.  Nonlinear  Rayleigh  Waves 

II.  Finite  Amplitude  Pulses  in  Liquids  with  Strong  Absorption 

III.  Reflection  of  a  Focused  Beam  from  a  Curved  Target 

IV.  Finite  Amplitude  Propagation  in  a  Sound  Velocity  Channel 

V.  Second  Harmonic  Generation  in  Pekeris  Waveguides 

VI.  Finite  Amplitude  Propagation  in  Multiple  Waveguide  Modes 

VII.  Cooperative  Radiation  and  Scattering  of  Sound  by  Bubbles 

VIII.  Water  Fountains  Produced  by  Intense  Standing  Waves  in  Air 

IX.  Spark-Soiurce  Lithotripter  Pulses 

Contributions  to  these  projects  were  marie  by  the  following  individuzds: 

Senior  Personnel 

•  M.  F.  Hamilton,  principzd  investigator 

•  Yu.  A.  Il’insky,  visiting  scientist 

•  E.  A.  Zabolotskaya,  visiting  scientist 
Graduate  Students 

•  M.  A.  Averkiou,  Ph.D.  student  in  Mechanical  Engineering 

•  C.  E.  Bruch,  M.A.  student  in  Physics 

•  Y.-S.  Lee,  Ph.D.  student  in  Mechaniced  Engineering 

•  D.  E.  Reckamp,  M.S.  student  in  Mechsiniced  Engineering 

•  D.  J.  Shull,  M.S.  student  in  Electrical  Engineering 

•  T.  W.  VanDoren,  Ph.D.  student  in  Mechanic2d  Engineering 
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Professor  Il’insky  ^^ld  Academicicin  Zabolotskaya,  on  leave  from  the  Physics  De¬ 
partment  of  Moscow  University  cmd  from  the  General  Physics  Institute  in  Moscow, 
respectively,  joined  our  acoustics  group  in  September  1991.  Professor  Il’insky  re¬ 
turned  to  Moscow  in  June  1992.  Their  meiin  source  of  fin^mci^d  support  was  the 
David  cmd  Lucile  Packard  Foundation  Fellowship  for  Science  emd  Engineering.  The 
graduate  students  received  partial  financial  support  from  several  funding  agencies, 
as  noted  in  the  relevcmt  sections. 

The  following  mcmuscripts  emd  abstracts,  which  describe  work  supported  at 
least  in  part  by  ONR,  have  been  published  (or  submitted  for  publication)  since 
1  August  1991. 

Refereed  Journeds 

•  M.  F.  Hamilton  and  E.  A.  Zabolotskaya,  “Nonhnear  propagation  of  sound  in 
a  liquid  layer  between  a  rigid  and  a  free  surface,”  J.  Acoust.  Soc.  Am.  90, 
1048-1055  (1991). 

•  M.  F.  Hamilton,  “Comparison  of  three  tremsient  solutions  for  the  ajcial  pres¬ 
sure  in  a  focused  sound  beam,”  J.  Acoust.  Soc.  Am.  92,  527-532  (1992). 

•  C.  M.  Daxvennes  emd  M.  F.  Hamilton,  “Additional  remarks  on  parametric 
reception  near  a  reflecting  pleme,”  J.  Acoust.  Soc.  Am.  (in  press). 

•  Yu.  A.  Il’insky,  emd  E.  A.  Zabolotskaya,  “Cooperative  radiation  emd  scatter¬ 
ing  of  acoustic  waves  by  gas  bubbles  in  liquids,”  J.  Acoust.  Soc.  Am.  (in 
press). 

•  M.  F.  Hamilton,  “A  tremsient  ajcial  solution  for  the  reflection  of  a  spherical 
wave  from  a  concave  ellipsoidal  mirror,”  J.  Acoust.  Soc.  Am.  (in  review). 

•  M.  F.  Hamilton,  Yu.  A.  Il’insky,  emd  E.  A.  Zabolotskaya,  “On  the  existence 
of  stationary  nonlinear  Rayleigh  waves,”  J.  Acoust.  Soc.  Am.  (in  review). 

Conference  Proceedings 

•  M.  F.  Hamilton,  “A  tremsient  solution  for  the  axial  pressure  field  of  a  spaurk- 
sovurce  lithotripter,”  to  appear  in  Proceedings  of  the  14th  International  Con¬ 
gress  on  Acoustics  (Beijing,  China,  September  1992). 

OrEil  Presentation  Abstracts 

•  M.  F.  Hamilton,  “An  ajcial  solution  for  the  reflection  of  a  spherical  wave  from 
a  concave  ellipsoidal  mirror,”  J.  Acoust.  Soc.  Am.  90,  2340  (1991). 

•  T.  W.  VcmDoren  emd  M.  F.  Hamilton,  “Water  fountains  produced  by  intense 
stemding  waves  in  Eiir,”  J.  Acoust.  Soc.  Am.  91,  2330  (1992). 
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•  Yu.  A.  Il’insky  and  E.  A.  Zabolotskaya,  “Cooperative  rzuiiation  of  acoustic 
waves  by  gas  bubbles  in  a  liquid,”  J.  Acoust.  Soc.  Am.  91,  2351  (1992). 

•  D.  E.  Reckamp,  E.  A.  Zabolotskaya,  and  M.  F.  Hamilton,  “Propagation  of 
finite  amplitude  sound  in  a  waveguide  with  a  parabolic  sound  velocity  profile,” 
J.  Acoust.  Soc.  Am.  91,  2352  (1992). 

•  M.  A.  Averkiou,  Y.-S.  Lee,  and  M.  F.  Hamilton,  “Propagation  of  pulsed 
finite  amplitude  sound  beams  in  a  liquid  with  strong  absorption,”  J.  Acoust. 
Soc.  Am.  91,  2455  (1992). 

•  M.  A.  Averkiou  and  M.  F.  Hamilton,  “Reflection  of  focused  sound  beams 
from  curved  surfaces,”  J.  Acoust.  Soc.  Am.  91,  2470  (1992). 

Theses 

•  D.  E.  Reckamp,  “Propagation  of  finite  amplitude  sound  in  a  waveguide  with 
a  parabolic  sound  velocity  profile,”  M.S.  Thesis,  The  University  of  Texais  at 
Austin  (May  1992). 


I.  Nonlinear  Rayleigh  Waves 


This  work  was  performed  by  Yu.  A.  Il’insky,  D.  J.  Shull,  and  E.  A.  Zabolotskaya. 
Additional  financieil  support  weis  provided  by  the  Packard  Foundation.  Comput¬ 
ing  resources  were  provided  by  The  University  of  Texeis  System  Center  for  High 
Performance  Computing. 

Rayleigh  waves  are  combinations  of  P  waves  (compressional  waves)  and  S  waves 
(shear  waves)  that  propagate  along  free  surfaces  of  elastic  solids.  Because  their 
energy  is  confined  within  a  layer  having  a  thickness  of  approximately  one  wave¬ 
length,  Rayleigh  waves  experience  relatively  small  spreading  losses  £is  compared 
with  unbounded  P  and  S  waves  (i.e.,  cylindrical  as  opposed  to  spherical  spread¬ 
ing).  Rayleigh  waves  therefore  transport  energy  from  earthquakes  and  underground 
explosions  over  large  distances  in  the  earth,  often  several  times  around.  Also,  the 
ease  with  which  Rayleigh  waves  cein  be  generated  and  detected  by  electrodes  on 
the  surfaces  of  piezoelectric  materieils,  and  the  faet  that  their  propagation  speed 
is  independent  of  frequency,  have  led  to  the  widespread  use  of  Rayleigh  waves  for 
performing  various  einaJog  signal  processing  functions  (e.g.,  delay  lines,  filters,  eind 
Fourier  transforms)  at  megahertz  frequencies. 

A  new  theoretical  model  for  the  propagation  of  nonlinear  Rayleigh  waves  has 
been  developed  recently  by  Zabolotskaya.*  The  model  is  based  on  Hamiltonizin 
formalism,  and  it  yields  a  set  of  coupled  equations  that  may  be  written  in  the 
dimensionless  form 


^  +  (^  +  ^n)K  =  2n^  i;  "i:  (1) 

an  \2K  J  m=n+l  m=l 
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where  Ki  is  a  complex  spectral  amplitude  in  a  Fourier  series  representation  of  the 
velocity  waveform,  R  is  distance,  An  are  frequency  dependent  attenuation  coeffi¬ 
cients,  and  Cmn  are  nonlinearity  coefficients  that  depend  on  the  individual  spectral 
interactions.  Equation  (1)  applies  to  plane  waves  for  m  =  0  and  to  cylindrical  waves 
for  m  =  1.  Three  investigations  based  on  Eq.  (1)  are  discussed  in  this  section. 

The  first  investigation  concentrates  on  the  properties  of  harmonic  generation 
in  nonlineeu:  Rayleigh  waves.  The  previous  malysis  by  Zabolotskaya*  focused  on 
distortion  of  the  time  waveform.  Quasilinear  solutions  have  now  been  derived  for 
second  harmonic  generation  in  plane  juid  cylindrical  Rayleigh  waves  that  propagate 
in  isotropic  solids  with  axbitreiry  attenuation  properties.  The  quasilinear  solutions 
were  used  to  derive  expressions  for  the  finite  zimplitude  losses  at  the  source  fre¬ 
quency  (following  the  very  successful  analysis  by  Merklinger  for  finite  amplitude 
soimd  in  thermoviscous  fluids^).  For  brevity,  we  present  here  only  the  approximate 
analytical  solutions  for  the  source  frequency  component  for  the  case  of  negligible 
sm2dl  signal  attenuation  (i.e.,  for  which  losses  occur  mainly  at  the  shocks): 


VI 


I  m  =  0 

^1  -I- 

_ _ 

v/l  +  32C?,i?o(v^- 


(plane) 


m  =  1 


(cylindrical) 


(2) 

(3) 


The  source  is  at  il  =  0  for  plane  waves  Jind  R  =  Rc  for  cylindrical  waves.  The 
above  expressions  are  found  to  be  in  good  agreement  with  the  numerical  solutions 
obtained  from  Eq.  (1),  even  well  into  the  shock  wave  region.  Note  that  the  source 
frequency  components  for  both  plane  and  cylindrical  waves  in  the  shock  wave  region 
(large  R)  axe  attenuated  at  the  same  rate  (^ls  i2“*).  Both  plane  eind  cylindrical 
Rayleigh  waves  axe  found  to  experience  saturation,  whereby  the  wave  amplitude  at 
a  fixed  field  point  no  longer  increases  with  source  level  for  source  levels  sufficiently 
high  that  shock  waves  axe  produced  near  the  source. 

In  Fig.  1  is  shown  the  distortion  and  eventual  acoustic  saturation  of  an  initially 
monofrequency,  planar  Rayleigh  wave  in  steel  as  a  function  of  source  level.  The 
spectra  and  waveforms  axe  normalized  by  their  values  at  the  source,  and  the  field 
point  is  at  a  fixed  distance  from  the  source.  Figiires  in  the  same  row  correspond  to 
the  same  normalized  source  level  It),  aJid  successive  rows  correspond  to  increasing 
source  level.  The  first  colunm  reveals  the  harmonic  generation  that  accompanies 
propagation  as  the  source  level  is  increased.  The  second  coliuim  reveals  the  dis¬ 
tortion  of  the  horizontal  velocity  waveform.  At  Lq  =  30  dB  the  waveform  exhibits 
the  sawtooth  profile  that  is  chaxacteristic  of  finite  amplitude  sound  in  fluids.  At 
Lq  =  40  dB,  however,  cusps  appear  at  the  beginning  and  end  of  the  shock  front. 
The  peak  amplitudes  at  these  cusps  generziUy  exceed  the  amplitude  of  the  source 
waveform.  Since  the  depth  penetration  of  a  Rayleigh  wave  is  inversely  proportional 
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Figtire  1:  Predicted  frequency  spectra  and  time  waveforms  associated  with  the 
nonlinear  distortion  and  acoustic  saturation  of  a  Rayleigh  wave,  as  a  function  of 
source  level. 
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to  frequency,  the  harmonic  generation  that  accompanies  waveform  distortion  tends 
to  localize  the  total  energy  in  the  wave  increasingly  closer  to  the  surface.*  When 
the  source  level  is  increased  another  10  dB,  to  Lq  =  50  dB,  saturation  has  oc¬ 
curred,  and  the  wave  amplitude  at  the  field  point  is  now  less  than  that  at  the 
source  because  of  increased  losses  at  the  shocks.  Further  increase  in  source  level 
merely  increases  these  losses,  but  the  profile  of  the  waveform  remains  the  same.  In 
the  third  column  axe  shown  the  vertical  velocity  waveforms  (which  axe  the  Hilbert 
transforms  of  the  horizontal  velocity  waveforms),  and  in  the  fourth  colunm  axe 
shown  the  vertical  displacement  waveforms  (which  axe  the  integrsds  of  the  verticzd 
velocity  waveforms). 

The  second  investigation  examines  whether  stationary  nonlinear  planar  Rayleigh 
waves  [i.e.,  for  which  dVn/dR  =  0  with  m  =  0  and  =  0  in  Eq.  (1)]  can  exist.^ 
Their  existence  has  been  predicted  by  Parker  and  Talbot,^"^  who  suggest  that  the 
nonlocal  character  of  nonlinear  Rayleigh  waves  provides  a  possible  mechanism  for 
the  existence  of  nondistorting  waveforms.  Parker  and  Talbot  used  numerical  solu¬ 
tions  of  a  model  equation  that  is  similar  in  form  to  Eq.  (1)  to  support  the  proposed 
existence  of  nonlinear  Rayleigh  waves.  However,  we  have  found  that  results  beised 
on  Eq.  (1)  with  dVn/dR  =  0  and  m  =  0  do  not  support  their  conclusions. 

Figure  2(a)  shows  a  typical  computed  “stationary”  spectrum  which  weis  limited 
to  iV  =  50  heirmonics.  The  corresponding  horizontal  velocity  waveform  is  shown 
in  Fig.  2(b)  (which  reveals  Gibbs  oscillations  due  to  the  abrupt  truncation  of  the 
frequency  spectrum).  If  N  can  be  chosen  sufficiently  large  that  no  further  increase 
in  N  introduces  a  change  in  either  the  spectrum  or  the  waveform,  then  one  may 
claim  that  a  stationary  waveform  has  been  identified  numerically  (this  was  the  pro¬ 
cedure  followed  by  Parker  and  Talbot).  However,  this  limit  can  never  be  achieved, 
as  was  proved  both  analytically  and  numerically  in  Ref.  3.  For  example,  doubling 
N  to  2N  will  produce  a  spectrum  with  exactly  the  same  spectrzd  envelope  as  in 
Fig.  2(a)  (when  renorm^llized  to  the  harmonic  scale  1  to  2N)  but  with  twice  the 
spectral  density.  In  other  words,  there  will  still  be  three  "spectral  nulls,”  but  they 
will  now  be  evenly  speiced  between  the  first  and  hundredth  harmonic,  rather  than 
the  first  and  fiftieth.  The  waveform  obviously  changes  in  this  ceise  (the  pulse-like 
shape  of  the  waveform  becomes  narrower  because  of  the  higher  frequency  content). 
In  conclusion,  a  stationzuy  wave  solution  cannot  be  obtjiined  on  the  basis  of  Eq.  (1). 

The  investigation  of  stationary  waves  is  completed  and  the  results  have  been 
submitted  for  publication.^  The  work  associated  with  Fig.  1  will  be  submitted  for 
publication,  and  it  is  currently  being  extended  to  diffracting  Rayleigh  wave  beams. 

The  third  investigation  produced  a  major  new  result.  In  particular,  a  method 
was  discovered  by  which  Eq.  (1)  can  be  rewritten  in  the  time  domain  £is  an  evolution 
equation  (analogous  to  the  Burgers  equation  for  finite  amplitude  sound  in  fluids). 
This  work  is  still  in  progress,  and  therefore  details  will  be  presented  in  the  Fifth 
Annual  Summary  Report  (i.e.,  next  year).  Here,  we  identify  only  the  main  features 
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Figure  2:  NumericaJ  solution  for  a  stationary  wave  with  three  spectrsd  nodes,  for 
m  =  An  =  Q  and  N  —  50:  (a)  frequency  spectrum  2ind  (b)  horizontal  velocity 
waveform. 


of  the  result.  The  nonlinearity  coefficient  matrix  may  be  written  as  the  summation^ 
C'mn  =  C'mn  +  ^mn  +  Each  of  the  three  terms  zdone  produces  nonlinear 

distortion  effects  which  lead  to  waveforms  having  the  saune  profiles  as  in  Fig.  1. 
The  nonlinear  evolution  equation  eissociated  with  the  third  term  is  the  easiest  to 
derive,  and  it  is  given  in  dimensionless  form  by 


dD  mD 


-A 


dD  dD* 


+  D 


d^D 

dT^ 


(4) 


where  D*  is  the  complex  conjugate  of  the  displacement  veiriable  D,  T  is  time,  A 
is  an  attenuation  coefficient,  and  7{  is  the  Hilbert  treinsform  operator.  The  re2d 
part  of  the  complex  function  D  is  proportion2d  to  the  horizontzd  component  of 
the  particle  displacement,  and  the  imaginsuy  part  is  proportionzd  to  the  verticed 
component.  A  computer  progreim  has  been  written  to  solve  Eq.  (4),  and  the  results 
are  indistinguishable  from  those  obtained  by  solving  Eq.  (1)  with  Cmn  =  C^n 
An  =  n^A\.  It  can  be  shown  that  when  the  Burgers  equation  is  rewritten  in  terms 
of  the  complex  displacement  variable  D,  the  result  is  identical  in  form  to  Eq.  (4) 
apsirt  from  the  absence  of  the  term  D{d^D* /dT^).  This  term  accounts  for  the 
nonlocal  nonhnearity  associated  with  Rayleigh  waves  (which  is  responsible  for  the 
cusping  near  the  shocks  in  the  horizontzd  velocity  waveforms  in  Fig.  1). 


II.  Finite  Amplitude  Pulses  in  Liquids  with  Strong  Ab¬ 
sorption 

The  numeric2d  work  for  this  project  was  performed  by  Y.-S.  Lee,  the  experimen- 
t2d  work  by  M.  A.  Averkiou.  Additional  financial  support  was  provided  by  the 
Packard  Foundation  and  the  National  Science  Foundation.  Computing  resources 
were  provided  by  The  University  of  Texas  System  Center  for  High  Performance 
Computing. 
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The  numerical  work  by  Lee  was  described  briefly  in  the  Third  Annual  Summary 
Report®  under  the  present  ONR  grrint,  and  subsequently  reported  in  the  Proceed¬ 
ings  of  Ultrasonics  International  91.®  The  following  abstract  appeared  at  the  May 
1992  Meeting  of  the  Acoustical  Society  of  America  in  Salt  Lake  City:^° 

Measurements  of  intense  acoustic  pulses  generated  by  directive  sources 
in  a  liquid  with  strong  absorption  are  compared  with  theoretical  predic¬ 
tions  based  on  a  time  domain  numerical  solution  of  the  KZK  nonlinear 
parabolic  wave  equation  [obtained  previously  by  Lee  and  Hamilton®]. 

The  computer  program  is  useful  for  describing  the  propagation  of  tone 
bursts  with  various  amplitude  emd  frequency  modulations.  Experiments 
were  performed  in  glycerin  with  narrowbEind  pulses  having  center  fre¬ 
quencies  of  several  megahertz.  The  pulses  experience  self-demodulation 
that  leads  to  farfield  waveforms  characterized  by  the  low  frequency 
envelope,  as  predicted  by  the  asymptotic  theory  of  Berktay.^^  Good 
agreement  between  theory  and  experiment  w£is  obtained  for  amplitude 
modulated  tone  bursts,  not  only  at  farfield  axial  positions  ^ls  achieved 
first  by  Moffett  et  al.,^^  but  also  off  axis,  throughout  the  tremsition  re¬ 
gion  emd  into  the  nearfield.  The  effects  of  frequency  modulation  are 
also  discussed. 

An  acceleration  under  this  ONR  gremt  provided  for  the  purchase  of  a  LeCroy 
9112  arbitrEiry  function  generator  (12  bits  dynamic  remge,  20  nsec/point  output 
rate).  This  high  precision  function  generator  has  made  possible  a  detailed  study  of 
nonlineeir  effects  in  pulsed  sound  beams.  An  example  of  experiment  emd  theory  is 
presented  in  Fig.  3.  The  pulse  has  a  center  frequency  of  3.5  MHz  emd  was  produced 
with  a  1.3  cm  diameter  source.  Glycerin  wzis  selected  8is  the  medium  because  of 
its  high  absorption  (approximately  6  dB/cm  at  3.5  MHz),  which  is  required  to 
investigate  self-demodulation  within  the  dimensions  of  the  ultr^lsonic  temk  facility 
in  the  MecheiniceJ  Engineering  Department.  The  decibel  levels  cited  next  to  each 
waveform  indicate  level  relative  to  that  at  the  source,  emd  the  dimensionless  dis¬ 
tance  a  =  zfzQ  is  in  terms  of  the  Rayleigh  distemce  Zq  at  3.5  MHz.  Results  from 
theory  aind  experiment  are  seen  to  be  in  outstanding  agreement,  which  attests  to 
the  precision  of  both  the  Mechanical  Engineering  teink  facility  aind  the  numerical 
code.  The  slight  asymmetry  of  the  final  experimented  waveform  is  attributed  to 
slight  asymmetry  in  the  response  of  the  source  tremsducer  (neimely,  the  ringing 
associated  with  the  high  Q),  which  was  not  taken  into  account  in  the  theoretical 
predictions.  The  results  in  Fig.  3  represent  the  first  instance  in  which  theory  and 
experiment  for  self-demodulation  have  been  compared  at  dist^mces  before  the  final 
demodulated  waveform  (in  the  bottom  row)  has  been  established. 
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III.  Reflection  of  a  Focused  Beam  from  a  Curved  Target 

This  work  was  performed  by  M.  A.  Averkiou,  who  received  additional  financial 
support  from  the  Packard  Foundation  and  the  National  Science  Foundation.  The 
work  was  reported  at  the  May  1992  Meeting  of  the  Acoustical  Society  of  America 
in  Sedt  Lake  City.^^  The  abstract  of  that  presentation  follows: 

The  reflection  of  focused  sound  beams  from  surfaces  with  sphericed 
curvature  is  investigated  theoretically  and  experimentadly.  Theoreti- 
ced  predictions  for  the  incident  and  reflected  beams  are  based  on  the 
p^lrabolic  wave  equation.  A  circul^lr  source  with  a  uniform  amplitude 
and  quadratic  phase  distribution  is  assumed.  Solutions  for  the  reflected 
beam  are  derived  for  both  pulsed  eind  continuous  sources.  The  experi¬ 
ments  were  performed  in  water  with  a  3.5  MHz  source  that  h^ls  a  nomi¬ 
nal  radius  of  2.5  cm  ^^ld  foced  length  of  15  cm.  Accurate  measurements 
of  the  incident  beam,  particularly  very  near  the  source,  were  used  to 
chairacterize  the  effective  radius  ^lnd  foced  length.  Reflection  from  both 
convex  and  concave  surfaces  weis  investigated.  The  targets  were  made 
of  nickel  with  radii  of  curvature  that  vary  from  5  cm  up  to  infinity  (pla- 
n£ir  tEirgets).  Measurements  of  the  reflected  beam  were  obtained  with 
a  needle  hydrophone  that  passed  through  a  small  hole  in  the  center  of 
the  source.  Agreement  between  theory  and  experiment  is  excellent,  and 
the  results  suggest  novel  ways  to  measure  surface  curvature. 

The  theory  is  based  on  a  solution  for  the  reflection  of  a  pulsed,  focused  sound 
beam  from  a  rigid  spherical  teirget.  For  simplicity,  we  present  here  only  the  solution 
for  the  case  in  which  the  ajds  of  the  be£im  pEisses  through  the  center  of  the  target: 

n  /rv,  -  /(^)  -  /(^  -  '^r)  .  . 

l-zld  +  2h-^{z-Zr){l-Zrld)  ^  ^ 

where  po  is  the  pressure  amplitude  at  the  source,  /(<)  is  the  source  waveform, 
r  =  t  —  zjcQ  is  retarded  time  in  terms  of  the  coordinate  z  along  the  axis  of  the 
beEun,  d  is  the  focal  length  of  the  source,  h  is  the  radius  of  ciirvature  of  the  target, 
z  =  Zr  is  where  the  axis  of  the  beam  intersects  the  target,  Eind 

_  fl  —  zld  +  2b~^(z  —  Zr)(l  —  Zr/d) 

2c^  \  Z  -f  26~^Zr(^  —  ^r) 

where  a  is  the  source  radius.  The  solution  is  written  in  terms  of  an  “unfolded 
geometry,”  in  which  the  beam  is  assumed  to  propagate  in  the  positive  z  direction 
following  reflection  from  the  target.  This  geometry  also  facilitates  the  presentation 
of  data.  For  example,  for  a  planar  target  (6  =  oo)  the  solution  for  pr  reduces  to  the 
freefleld  solution  for  a  focused  beam  in  the  parabolic  approximation.  Modification 
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of  Eq.  (5)  to  account  for  a  geometry  in  which  the  center  of  the  target  does  not 
coincide  with  the  beam  axis  is  straightforward. 

Shown  in  Fig.  4  are  comparisons  of  theory  with  experiment  for  the  amplitude  of 
a  sinusoidal  source  waveform  that  is  reflected  from  either  a  convex  target  (6  =  +ci) 
or  a  concave  target  (6  =  —d).  The  experiment  was  performed  in  water,  and  the 
reflected  pressure  was  measured  in  the  center  of  the  source,  which  corresponds  to 
the  location  z  =  2zr  (i.e.,  twice  the  round  trip  distance  from  the  source  to  the 
target).  The  receiver  consisted  of  a  small  needle  hydrophone  inserted  in  the  center 
of  the  source  (the  source  radius  was  a  =  2.5  cm,  and  its  focal  length  was  d  = 
15.2  cm).  The  vertical  ajces  in  Fig.  4  correspond  to  pressure  amplitude  (majcimized 
by  its  maximum  value),  the  horizontal  axes  to  the  distance  between  the  source 
and  the  target.  If  the  target  were  planar,  the  received  pressure  amplitude  would 
be  majcimized  with  a  source/target  separation  given  by  Zrjd  =  0.5,  i.e.,  when  the 
target  is  one  half  focal  length  away  from  the  source  (because  the  received  signal  h^ls 
traveled  a  roimd-trip  distance  d).  The  arrows  in  the  figure  indicate  this  position 
and  provide  a  point  of  reference  for  comparison  with  results  for  curved  targets.  For 
the  concave  target  (6  =  — d),  the  maximum  is  measured  at  Zr/d  <  0.5  (because 
the  target  enhances  the  focusing  effect  and  thus  causes  the  beam  to  focus  within 
a  shorter  distance,  i.e.,  the  net  focal  length  decreases).  For  the  convex  target 
(6  =  +d)  the  maximum  is  measured  at  Zr/d  >  0.5  (because  the  target  defocuses 
the  beam  and  thus  increases  the  net  focal  length).  Theory  and  experiment  are 
seen  to  be  in  very  good  agreement.  The  oscillations  in  the  response  curves  are 
due  to  diffraction,  specifically,  the  interference  of  the  center  wave  and  edge  wave 
radiated  by  the  source.  These  contributions  are  represented  by  the  two  terms  in 
the  numerator  of  Eq.  (5). 


Figure  4:  Q)mparison  of  linear  theory  and  experiment  for  the  reflection  of  a  focused 
sound  beeun  from  either  a  concave  target  (6  =  —d)  or  a  convex  tairget  (d  =  +&). 
The  point  receiver  is  mounted  in  the  center  of  the  source  trainsducer. 

The  experiments  reported  here  will  be  extended  to  the  more  general  C2ise  of 
transmission  eind  reflection  at  a  curved  fluid-fluid  interface,  for  both  finite  ampli- 
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tude  as  well  as  for  small  signal  sound  beams. 


IV.  Finite  Amplitude  Propagation  in  a  Sound  Velocity  Chan¬ 
nel 


The  personnel  working  on  this  project  were  Ensign  D.  E.  Reckamp,  USN,  and  E.  A. 
Zabolotskaya.  The  main  ^in^lnciaI  support  for  Reckamp  was  provided  by  an  Applied 
Research  Laboratories  Naval  Academy  Scholarship.  Reckcimp  reported  this  work 
in  his  M.S.  thesis. He  received  an  M.S.  degree  in  Mechzmiced  Engineering  in  May 
1992.  Apart  pubhcation  of  the  results  in  a  journal  article,  this  project  is  completed. 
Additional  fin2mcial  support  weis  provided  by  the  Packard  Foundation. 

Background  on  the  project  is  discussed  in  the  Third  Annual  Summary  Report* 
under  the  present  ONR  grant.  The  results  were  presented  at  the  May  1992  Meet¬ 
ing  of  the  Acoustical  Society  of  America  in  Salt  Leike  City.  The  following  is  the 
abstract**  that  accompanied  the  presentation; 


The  propagation  of  finite  amplitude  sound  in  a  waveguide  with  a  para¬ 
bolic  sound  velocity  profile  is  investigated  theoretically.  It  is  sissumed 
that  the  prim^lry  wave  propagates  in  a  single  mode  at  a  frequency  that  is 
large  compared  with  the  cutoff  frequency.  The  acoustic  energy  is  there¬ 
fore  concentrated  near  the  axis  of  the  sound  channel.  Both  the  primary 
wave  and  the  nonlinearly  generated  second  harmonic  component  have 
mode  shapes  that  are  described  by  Gauss-Hermite  eigenfunctions.  For  a 
primary  wave  in  mode  m,  the  second  harmonic  component  is  generated 
in  8l11  even  order  modes  up  to  2m,  with  most  of  the  energy  contmned  in 
mode  2m.  A  nonlineeu:  Schrodinger  equation  is  derived  for  the  envelope 
of  a  neurrowband  pulse  that  propagates  in  a  single  mode.  Analytical 
expressions  eure  derived  for  the  coefficients  in  the  Schrodinger  equation, 
zmd  it  is  found  that  only  dark  envelope  solitons  can  propagate  without 
distortion  in  the  waveguide. 


The  mode  shapes  of  the  second  harmonic  component  were  presented  in  the  Third 
AnnueJ  Snmma.ry  Report.*  Below,  we  present  the  general  results. 

The  governing  wave  equation  is  the  following  modified  Westervelt  equation: 


vV- 


dt^ 


J^ 

PoCq 


(7) 


where  z  is  depth  in  the  water  colunm  and  fi  chEU’acterizes  the  curvature  of  a  peu’a- 
bolic  sound  velocity  profile  about  the  waveguide  eods  at  r  =  0  (see  Fig.  5).  This 
sound  speed  profile  is  frequently  used  to  model  the  SOFAR  channel.  The  limita¬ 
tions  of  Eq.  (7)  are  that  the  acoustic  wavelengths  miist  be  short  compared  with 
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Figure  5:  Coordinate 


system  for  the  parabolic  sound  velocity  channel. 


the  characteristic  dimension  of  the  channel  (i.e.,  the  effective  width  of  the  waveg¬ 
uide),  and  that  the  sound  must  be  trapped  near  the  waveguide  eixis  (which  means 
the  source  frequency  must  be  severed  times  the  cutoff  frequency).  If  it  is  further 
assumed  that  the  primeiry  sound  field  consists  of  a  single  frequency  u  and  propa¬ 
gates  in  a  single  mode  m,  the  following  quasilinear  solutions  axe  obtained  for  the 
fundeimental  (pi)  and  second  harmonic  (pj)  pressure  components: 

p.  =  (8) 

where  is  the  wavenumber  of  the  fundamentad  component,  2/?^  is  the  axisd 
wavenumber  of  a  freely  propagating  second  harmonic  component,  Eind  is  the 
Hermite  polynomial  of  order  m.  The  expression  for  pi  is  well  known,*®  whereas 
the  expression  for  pj  is  new.  Because  the  fundamentad  and  second  harmonic  com¬ 
ponents  propagate  at  different  speeds  in  the  waveguide  {km  ^  «m).  the  amplitude 
of  the  second  harmonic  component  experiences  spatizd  beating.  The  spatial  peri¬ 
odicity  of  the  beating  is  7r/(Km  ~  ^m)  —  47r/pci).  The  expression  for  pj  is  a  single 
mode  approximation.  The  full  expression  for  pj  includes  m  -|-  1  modes,  i.e.,  £dl  even 
modes  from  0  through  2m.  However,  it  W8is  discovered  that  most  of  the  energy  is 
contained  in  the  highest  mode,  2m,  which  provides  justification  for  the  single  mode 
approximation. 
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The  propagation  of  narrowband  pulses  was  also  considered.  The  expression  for 
a  narrowband  pulse  that  propagates  in  a  single  mode  may  be  written 


(10) 


where  ujq  is  the  center  frequency  of  the  pulse,  a  is  the  pulse  envelope  function, 
T  =  t  —  xfCg  IS  a  reteirded  time,  and  Cg  the  group  velocity  of  the  pulse.  A  nonlineair 
Schrodinger  equation  (NLS)  can  be  derived  for  waveguides:*^ 


’dx  2 


+  A|a|^a  =  0 


(11) 


where  k'^  =  (Pk^/duQ.  The  NLS  predicts  the  existence  of  envelope  solitons,  bright 
solitons  if  A./k'^  >  0  or  dark  solitons  if  A/ A:"  <  0.  The  case  A/A:^  <  0  also  indicates 
that  localized  energy  peickets  (i.e.,  which  vanish  as  r  — >  ±oo)  will  disperse  as  they 
propagate.  An  analytical  expression  was  derived  for  the  ratio  A/ A:"  in  the  single 
mode  approximation  for  the  second  heirmonic  component:*^ 


where  u!m  is  the  cutoff  frequency  for  the  pulse  (ti;o  >  u>m  is  eissumed).  This  ex¬ 
pression  is  always  negative,  which  means  that  bright  solitons  cannot  exist  in  this 
waveguide,  and  pulses  will  disperse  £is  they  propagate. 


V.  Second  Harmonic  Generation  in  Pekeris  Waveguides 

This  work  was  performed  by  C.  E.  Bruch,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya, 
who  each  received  additional  financial  support  from  the  Packard  Foimdation.  Fur¬ 
ther  background  information  on  this  project  may  be  found  in  the  Third  Annual 
Summary  Report.® 

The  Pekeris**  waveguide  model  consists  of  two  liquid  layers  having  different 
soimd  speeds.  The  upper  layer  (0  <  z  <  hi)  is  boimded  above  by  a  free  surface 
at  2  =  0  and  below  (z  >  hi)  by  a  second,  penetrable  layer.  The  lower  layer  is 
usuailly  assinned  to  have  infinite  depth  (h2  =  oo).  When  the  sound  speed  in  the 
lower  layer  exceeds  that  in  the  upper  layer  (c2  >  Ci,  where  the  index  1  refers  to  the 
upper  layer  emd  the  index  2  to  the  lower  layer),  guided  wave  propagation  can  exist 
in  the  upper  layer.  This  is  the  situation  for  sea  water  over  sediment.  Although  the 
Pekeris  waveguide  model  has  been  widely  used  to  study  the  zicoustics  of  shallow 
ocean  channels  in  the  small  signal  approximation,*®  no  analyses  have  been  reported 
on  finite  amplitude  effects.  Here  we  present  preliminary  results  for  second  harmonic 
generation  in  a  Pekeris  waveguide. 

In  the  course  of  the  investigation,  it  proved  useful  to  consider  a  slightly  modified 
version  of  the  Pekeris  waveguide  model.  In  particular,  a  lower  layer  of  finite,  rather 
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than  infinite,  depth  weis  assumed  (/ii  <  /12  <  00).  The  reason  for  this  assumption 
is  that  the  sound  field  in  Ein  infinite  domain  (/12  =  00)  must  be  described  with  a 
continuous  spectrum  of  modes  (plus  a  discrete  spectrum  that  corresponds  to  the 
upper  layer),  whereas  a  field  in  a  finite  domain  consists  of  only  a  discrete  spectrum 
of  modes.  The  continuous  spectrum  is  frequently  ignored,  because  it  contributes 
significEintly  to  the  acoustic  field  only  very  neax  the  source.'®  However,  Einalysis 
of  second  hEirmonic  generation  requires  caxeful  consideration  of  the  field  necir  the 
source,  in  particulax,  to  properly  match  the  source  condition.  As  a  practical  matter, 
matching  the  source  condition  EuialyticEilly  is  more  easily  achieved  when  the  modes 
form  a  discrete  rather  than  a  continuous  spectrum.  For  the  case  of  trapped  wave 
propagation  in  the  upper  layer,  in  which  acoustic  penetration  in  the  lower  layer 
is  restricted  to  within  approximately  one  wavelength  beneath  the  interface,  it  was 
found  that  the  thickness  of  the  lower  layer  need  only  be  on  the  order  of  the  thickness 
of  the  upper  layer  (/12  ~  2/ii)  to  simulate  2Ln  ideal  Pekeris  waveguide  (/12  =  00). 


Mode  Shape 
(m=1,n=14.5) 


z/h. 


Mode  Shape 
(m=2,  n=32.6) 


Mode  Shape 
(m=3,  n=50.7) 


Second  Harmonic 

Particular  Solution; -  Unear  Solution: 


Figure  6:  Predicted  fundamental  and  second  haxmonic  mode  shapes  (dashed  and 
solid  lines,  respectively)  in  a  Pekeris  waveguide  formed  by  sea  water  over  sediment. 

In  Fig.  6  Eire  shown  mode  shapes  corresponding  to  second  haxmonic  generation 
(solid  lines)  by  a  primaxy  wave  (deished  lines)  that  propagates  in  mode  m  =  1, 
2,  or  3  in  a  layer  of  sea  water  over  sediment  (ci/c2  =  0.86,  pi/p2  =  0.5).  The 
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bottom  of  the  lower  layer  (/12  =  2/ii)  is  Eissumed  to  be  perfectly  rigid.  Only  the 
forced  component  (the  particulzir  solution)  of  the  second  heirmonic  wave  is  shown 
in  the  figure.  The  dimensionless  source  frequency  =  uc\  / /ij  increaises  with  mode 
number  in  such  a  way  eis  to  maintain  a  constzmt  zmgle  of  incidence  between  the 
water-sediment  interface  zmd  the  propagation  directions  of  the  pl^me  waves  that 
constitute  the  mode  at  the  source  frequency.  This  zmgle  is  ^  =  10°.  The  critic^J 
zmgle  for  this  case  is  =  30°,  which  means  that  totEd  intem^J  reflection  occurs  in 
the  upper  layer,  zmd  the  primary  wave  is  trapped.  The  mmn  observation  is  that 
the  penetration  depth  of  the  nonlinearly  generated  second  heirmonic  component 
is  slightly  less  thzm  at  the  source  frequency.  Otherwise,  the  properties  of  second 
harmonic  generation  due  to  a  primary  wave  trapped  in  the  upper  layer  of  a  Pekeris 
waveguide  Eire  similar  to  those  for  the  case  in  which  the  sound  in  the  upper  layer 
is  trapped  by  a  rigid  surface  z  =  h\  For  exeimple,  second  heirmonic  generation 
is  slightly  more  efficient  near  the  free  surface  than  near  the  fluid-fluid  interface. 

VI.  Finite  Amplitude  Propagation  in  Multiple  Waveguide 
Modes 

This  work  was  performed  by  T.  W.  YzuiDoren,  who  received  additionzd  financiEil 
support  from  a  Rockwell  Fellowship  Euid  the  Packard  Foundation.  Preliminzury 
work  on  the  project  was  discussed  in  the  Third  Annual  Summary  Report.® 

The  waveguide  investigations  described  in  the  previous  two  sections  of  this 
report  (Secs.  IV  Euid  V)  Eire  based  on  the  Eissumption  that  the  primary  wave  propa¬ 
gates  in  only  one  mode,  which  results  in  a  tremendous  simplification  of  the  axiEdysis. 
Except  in  csurefully  designed  experiments,  however,  this  is  usuEdly  not  the  CEise,  Euid 
the  primary  wave  propagates  in  a  vEiriety  of  modes.  The  purpose  of  the  present 
investigation  is  to  consider  how  the  propagation  of  a  primary  wave  in  more  thEin 
one  mode  affects  hEirmonic  generation  Emd  shock  formation.  The  EmEdysis  proceeds 
as  follows.  First,  a  quEisilincEir  ariEilyticEil  solution  of  a  Westervelt  nonlinesu:  wave 
equation  is  derived  for  second  hEirmonic  generation  in  a  rectEuigulsu:  duct  with  rigid 
walls,  with  the  primary  wave  propagating  in  Ein  Eurbitrsuy  number  of  modes.  Sec¬ 
ond,  the  corresponding  quEisilineEir  solution  is  derived  on  the  basis  of  the  pEirabolic 
(KZK)  approximation  of  the  nonlinesu:  wave  equation.  Third,  if  these  solutions 
(Westervelt  Emd  KZK)  compEire  well  with  each  other  Emd  with  experiment,  then 
existing  computer  programs  for  solving  the  fully  nonlinesu:  KZK  equation  will  be 
modified  for  application  to  guided  wave  propagation.^’^^  The  numericEil  solutions 
cEm  be  used  to  predict  third  Emd  higher  hEirmonic  generation,  as  well  sis  shock  for¬ 
mation.  Moreover,  pulsed  soimd  in  waveguides  can  then  be  explored.  Although  a 
numericEil  solution  of  a  KZK-type  equation  for  guided  finite  amplitude  sound  has 
Eilready  appesured  in  the  literature,^^  we  feel  that  a  careful  analysis,  bEised  on  both 
anEilyticEil  solutions  and  experiment,  is  needed  to  justify  this  approach. 
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Figure  7:  Comparison  of  quasilinear  theory  for  second  harmonic  generation  (solid 
lines)  with  experiment  for  a  33  kHz  piston  source  in  an  air-filled  rectangular  waveg¬ 
uide. 
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Quasilinear  analytical  solutions  for  second  harmonic  generation  by  a  primary 
wave  in  multiple  waveguide  modes  were  derived  from  both  the  Westervelt  and  KZK 
equations.  The  solutions  axe  too  cumbersome  to  present  here.  An  experiment  was 
performed  in  which  the  fundamentEd,  second  and  third  harmonic  components  were 
measured  in  am  air-filled  rectangular  duct.  The  transverse  dimension  of  the  duct  is 
6.6  cm,  and  it  was  excited  with  a  centered  3.3  cm,  33  kHz  source  that  (in  theory) 
produced  sound  in  four  propagating  modes.  The  second  harmonic  should  then 
be  generated  in  only  two  propagating  modes  for  this  geometry.  The  results,  both 
theory  (based  on  the  Westervelt  equation  with  ad  hoc  attenuation  coefficients) 
and  experiment,  are  shown  in  Fig.  7.  The  agreement  for  the  fundamental  and 
second  harmonic  components,  in  terms  of  general  properties,  is  very  reasonable. 
No  analytical  solution  for  the  third  heirmonic  component  was  derived.  Theoretical 
description  of  this  component  will  be  based  on  a  numerical  solution  of  the  KZK 
equation. 

VTI.  Cooperative  Radiation  and  Scattering  of  Sound  by 
Bubbles 

Work  on  this  project  was  performed  by  Il’insky  and  Zabolotskaya,  who  received 
additional  financial  support  from  the  Packard  Foimdation.  This  research  is  now 
completed,  auid  a  journal  article  has  been  accepted  for  publication  in  the  Journal 
of  the  Acoustical  Society  of  America. The  results  have  also  been  presented  at  the 
May  1992  Meeting  of  the  Acoustical  Society  of  America  in  SEilt  Lake  City.*'*  The 
following  is  the  abstract  from  that  meeting: 

The  process  of  cooperative  radiation  of  acoustic  waves  by  radial  bubble 
oscillations  is  investigated  theoretically.  In  addition  to  instemtemeous 
bubble  interaction,*®  we  have  also  teiken  into  account  the  interaction 
due  to  acoustic  radiation.  Numerical  results  for  the  acoustic  intensity 
produced  by  10  bubbles  that  pulsate  in  water  with  random  initial  phases 
are  presented.  Mutual  bubble  interaction  causes  synchronization  of 
the  bubble  pulsation  phases,  which  gives  rise  to  collective  radiation. 

The  phenomenon  is  not  as  strong  as  in  optics.  In  order  to  ensure 
superradiance  by  10  bubbles,  the  radial  oscillations  of  the  bubbles  must 
have  amplitudes  of  the  order  of  their  initial  radii. 

The  details  of  this  work  are  summarized  by  the  following  (modified)  excerpts 
from  the  forthcoming  joumrJ  article.*^ 

It  is  known  that  distributions  of  gas  bubbles  produce  dramatic  changes  in  the 
properties  of  a  liquid.  In  particular,  bubbles  introduce  strong  nonlinearity  and 
dispersion.  Bubble  oscillations  result  in  the  radiation  of  acoustic  waves,  eind  their 
contribution  to  oceam  noise  auid  sound  scattering  is  substamtiaJ.*® 
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The  aforementioned  phenomena  are  usually  analyzed  with  the  Rayleigh  equa¬ 
tion  for  single  bubble  pulsations.  Cooperative  effects  are  not  taken  into  accoimt, 
apart  from  the  bubble  contribution  to  the  coherent  part  of  the  acoustic  field  that 
arises  dining  sound  propagation  through  a  quasi-homogeneous  mixture  of  a  liquid 
with  bubbles.  This  is  the  usual  model  of  a  two-phase  homogeneous  medium.  The 
present  work  constitutes  a  theoretical  analysis  of  the  collective  radiation  and  scat¬ 
tering  of  sound  by  bubbles.  A  new  feature  of  the  analysis  is  the  inclusion  of  a  term 
that  takes  into  ^lccount  the  bubble  interaction  through  the  radiation  field. 

The  following  coupled  dimensionless  equations  were  derived  to  describe  the 
volume  oscillations  of  a  collection  of  bubbles  with  a  distribution  of  radii: 

dV 

(i3) 

m 

where  is  the  (dimensionless)  volinne  of  the  nth  bubble,  T  is  time,  accounts  for 
the  frequency  shift  connected  with  instantaneous  bubble  inter^lction  through  the 
radiation  field  and  the  deviation  of  the  bubble  oscillation  frequency  from  the  mean 
frequency  that  characterizes  the  given  bubble  distribution,  r„  is  a  coefficient  that 
accounts  for  both  quadratic  and  cubic  nonlinearity,  and  the  sinnmation  represents 
bubble  interaction  through  the  radiation  field. 


Figure  8:  Intensity  of  sound  radiated  by  ten  oscillating  bubbles  as  a  function  of 
time,  based  on  (a)  lineeir  and  (b)  nonlinear  theory.  The  standard  deviation  of  the 
bubble  radius  distribution  is  a. 


The  results  of  integrating  Eqs.  (13)  are  shown  in  Fig.  8.  A  field  of  10  bubbles  was 
assumed,  the  initial  pulsation  amplitudes  of  which  were  assumed  to  be  equal.  The 
initial  phase  of  each  bubble  is  Eirbitrary,  and  the  phase  distribution  (corresponding 
to  values  of  S„)  was  assigned  randomly  according  to  the  standard  deviations  a  = 
0, 2, 4, 6.  The  vertical  axes  in  Fig.  8  represent  I  =  1 which  is  proportional 
to  the  radiation  intensity.  Figure  8(a)  pertzdns  to  a  linear  (r„  =  0)  solution  of 
Eqs.  (13),  Fig.  8(b)  to  a  nonlineEu:  (r„  =  10)  solution.  Note  that  in  the  linear  case. 
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Fig.  8(a),  the  intensity  does  not  decay  monotonically  with  time  when  there  is  a 
distribution  of  bubble  sizes.  Locedized  intensity  maxima  occur  when  the  bubbles 
momentarily  oscillate  in  phase  as  a  result  of  their  different  resonance  frequencies. 

In  the  nonlinear  system,  Fig.  8(b),  the  locadized  intensity  maxima  are  an  order  of 
magnitude  greater  them  in  the  linear  ceise.  This  nonlinear  acoustical  phenomenon  is 
the  anedog  of  superradiance  in  nonlinear  optics.  To  achieve  nonlinearity  designated 
by  r„  =  10,  however,  the  bubbles  must  pulsate  with  volume  amplitudes  that  axe 
of  the  order  of  their  initied  volumes.  It  is  therefore  unlikely  that  a  superradiance 
pulse  will  be  observed  in  acoustics  via  bubble  interations  in  liquids. 

The  main  conclusion  of  this  work  is  that  cooperative  effects  should  be  taken  into 
account  when  there  axe  many  bubbles  in  the  liquid  and  they  axe  sufficiently  close  to 
each  other  to  interact.  Nonlinear  effects  cause  the  bubble  oscillation  frequencies  to 
depend  on  pidsation  amplitudes,  which  leads  to  phase  matching  and  a  phenomenon 
emalogous  to  superradiance  in  optics. 

VIII.  Water  Fountains  Produced  by  Intense  Standing  Waves 
in  Air 

This  work  was  primarily  experimented  and  was  performed  by  T.  W.  VanDoren. 
VanDoren  received  additioned  finemcied  support  through  a  Rockwell  Fellowship  and 
from  the  Peickard  Foundation.  The  following  abstract^  from  an  oral  presentation 
at  the  May  1992  Meeting  of  the  Acoustical  Society  of  America  in  Sedt  Lake  City 
svunmarizes  the  motivation  for  this  work. 

An  eicoustics  demonstration  found  in  severed  science  museiuns  involves 
a  horizonted  stemding  wave  tube  which  contains  a  shedlow  layer  of  liq¬ 
uid.  When  the  liquid  is  water  and  the  sound  pressvure  level  of  a  stand¬ 
ing  wave  in  the  air  above  the  water  exceeds  approximately  160  dB  (re 
20  /iPa)  at  frequencies  below  1  kHz,  vigorous  fountains  occirr  at  the 
nodes  in  the  pressure  field.  The  introduction  of  a  surfactemt  in  the 
water  causes  clusters  of  fountains  to  appear  near  the  nodal  planes  at 
somewhat  lower  sound  pressure  levels.  Standaxd  quasilineax  theory  pro¬ 
vides  a  reasonable  description  of  the  acoustic  streaming  and  dc  pressure 
in  the  air  prior  to  the  appearance  of  fountains.  The  foimtains  seem  to 
result  from  the  ejection  of  droplets  by  inner  streaming  vortices  that  axe 
formed  above  the  eur-water  interface.  The  water  fountains  and  stream¬ 
ing  patterns  axe  discussed  in  relation  to  Andrade’s  observations^^  of 
solid  paxticulate  motion  produced  by  intense  standing  waves  in  air. 

The  presentation  was  accompanied  by  a  videotaped  visualization  of  the  acoustic 
streeiming  in  the  standing  wave  tube.  A  half  cylinder  on  one  wjdl  of  the  tube, 
rather  than  the  small  modulations  produced  by  the  surface  of  the  water,  was  used 
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to  enhance  the  streaming.  The  vizualization  w^ls  performed  with  a  laser  and  smoke 
particles.  However,  the  inner  streaming  vortices  that  are  believed  to  be  responsible 
for  the  water  fountciins  could  not  be  vizualized  with  this  technique,  and  only  the 
outer  streaming  vortices  were  readily  viewed.  L^lser  doppler  velocimetry  was  used 
to  confirm  the  existence  of  the  inner  streaming  layer.  Whereas  acoustic  streaming 
patterns  of  this  type  (i.e.,  inner  £ind  outer  circulations  neeir  boundaries)  axe  well 
known,^*  ^ls  axe  water  fountains  produced  by  standing  waves  in  air  (they  may  be 
found  in  several  science  museums),  the  suggestion  that  the  fountains  axe  driven  by 
streeiming  has  not,  to  our  knowledge,  been  discussed  previously  in  the  literature. 


IX.  Spark-Source  Lithotripter  Pulses 

This  work  was  performed  mainly  during  the  1990-91  reporting  period,  as  described 
in  the  Third  Annual  Summary  Report®  under  the  present  ONR  grant,  but  the 
comparisons  of  theory  and  experiment  shown  below  axe  new.  The  project  is  now 
completed,  and  an  article  heis  been  submitted  for  publication  in  the  Journal  of  the 
Acoustical  Society  of  America.^®  This  work  was  also  presented  at  the  November 
1991  Meeting  of  the  Acoustical  Society  of  America  in  Houston,®®  and  it  is  scheduled 
for  presentation  at,  and  publication  in  the  proceedings  of,  the  14th  International 
Congress  on  Acoustics  to  be  held  in  Beijing  in  September  1992.®^  Additional  finan¬ 
cial  support  was  provided  by  the  Packard  Foundation. 

The  results  may  be  summarized  with  the  following  abstract  from  the  submitted 
journal  article:®® 

A  transient  solution  is  derived  for  the  refiected  pressure  field  along  the 
axis  of  symmetry  of  a  concave,  eixisymmetric,  ellipsoided  mirror.  The 
incident  field  is  a  spherical  wave  produced  at  the  near  focus  of  the 
mirror.  Short  wavelengths  in  comparison  with  the  minimum  radius  of 
curvature  of  the  mirror,  and  lossless,  small  signed  propagation  in  a  ho¬ 
mogeneous  fiuid  axe  eissumed.  The  pressure  at  the  feir  focus  is  given 
by  a  simple  anedytical  solution,  and  elsewhere  along  the  eixis  by  a  con¬ 
volution  integred.  The  results  constitute  a  generalization  of  the  theory 
developed  for  spherical  mirrors  by  Comet  and  Blackstock.®®  Applica¬ 
tion  to  waveforms  produced  by  spark-source  lithotripters  is  discussed, 
and  comp£irisons  eire  made  with  measurements  reported  by  Muller.®® 

The  axiad  solution  for  the  reflected  pressure  pa  can  be  written  in  the  form 

^  ^  r  H^(z,  (')/((  -  (')  dt'  (14) 

Po  a  •'<1 

where  po  is  the  pressure  aimplitude  and  f(t)  the  time  dependence  of  the  source 
waveform.  The  first  term  is  associated  with  the  center  wave,  the  second  term 
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with  the  edge  wave,  £ind  the  third  term  with  the  wake.  The  corresponding  impulse 
response  functions  He,  Hg  and  H,^  are  defined  in  Refs.  29  £ind  31. 

For  the  C£ilculations  in  Fig.  9  we  considered  a  half  ellipse  with  eccentricity 
e  =  0.7  (md  therefore  a  far  focus  at  z /a  =  0.7).  The  source  waveform  f[t)  is  shown 
in  Fig.  9(a).  The  paxauneters  approximate  those  of  m  experiment  performed  in 
water  by  Muller whose  measured  waveforms  (reproduced  from  Fig.  4  of  Ref.  33) 
at  three  axi8d  locations  {z/a  =  0.47,  0.75,  amd  0.93)  appear  in  the  left  colunm 
of  Fig.  9,  opposite  the  corresponding  theoretical  predictions  based  on  Eq.  (14). 
Muller  used  a  Dornier  XLl  lithotripter,  smd  the  duration  of  his  incident  pulse  was 
approximately  T  =  4  /isec.  The  pairs  of  vertic8il  dashed  lines  in  Figs.  9(c)-9(e) 
identify  the  approximate  beginnings  smd  ends  of  the  oscilloscope  traces  in  the  left 
colunm.  The  label  C  in  the  figures  identifies  the  beginning  of  the  center  wave,  E 
identifies  the  beginning  of  the  edge  wave,  amd  W  identifies  the  most  pronounced 
contribution  due  to  the  wake. 

At  z/a  =  0,  the  center  wave  amd  edge  wave  axe  well  sep^u^ated,  amd  the  wake 
produces  a  slight  negative  pressure  immediately  following  the  center  wave.  At 
z/a  —  0.47  the  center  wave  md  edge  wave  axe  still  clearly  resolved.  Comparison 
of  Fig.  9(c)  with  the  oscilloscope  trace  to  the  left  indicates  that  the  small  negative 
pressure  at  the  end  of  the  measured  waveform  is  evidently  due  to  the  wake  rather 
than  the  edge  wave.  The  oscilloscope  trace  appears  to  end  prior  to  the  sirrivaJ  of  the 
edge  wave.  At  x/a  =  0.75,  just  beyond  the  f^u^  focus,  the  center  wave  £ind  edge  wave 
overlap,  £ind  the  wake  produces  a  l^u^ge  positive  pressure.  However,  the  predicted 
waveform  in  Fig.  9(d)  appears  backw^u^ds  in  comparison  with  the  corresponding 
measured  waveform.  This  reversal  is  probably  due  to  nonlinear  effects  introduced 
by  the  l^u^ge  peak  pressure,  which  is  approximately  800  bar  (80  MPa).  Nonlinearity 
would  cause  point  W  in  Fig.  9(d)  to  catch  up  with  point  E  eind  thus  produce 
a  waveform  more  like  that  which  was  measured.  Note  also  that  the  amplitude 
ratio  of  point  W  in  Fig.  9(d)  to  point  C  in  Fig.  9(c)  matches  the  ratio  of  the 
corresponding  measured  pe^lk  pressures.  At  z/a  =  0.93,  both  the  measured  aind 
predicted  waveforms  possess  a  horizontal  plateau  that  follows  the  arrival  of  the 
edge  wave  (E),  leading  to  a  spike  produced  by  the  wEike  (W)  in  the  middle  of  the 
waveform.  Again,  nonlineeu:  effects  would  cause  point  W  in  the  predicted  waveform 
to  advance  in  time  relative  to  point  E.  Although  the  predicted  eimplitude  ratio  of 
point  W  to  point  E  matches  the  measured  ratio,  the  predicted  negative  pressure 
due  to  the  center  wave  is  much  greater  thain  was  measured.  At  x/a  =  2,  the  edge 
wave  has  sepzirated  from  the  center  wave,  and  the  effect  of  the  wake  is  reduced. 
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Figure  9:  Comparison  of  present  theory  (right  coliunn)  and  experiment  for 
lithotripsy  pulses  measured  in  water  by  Miiller.^ 
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